Relaxation Mechanisms in Glassy Dynamics: the Arrhenius and Fragile Regimes 
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Generic glass formers exhibit at least two characteristic changes in their relaxation behavior, first 
to an Arrhenius-type relaxation at some characteristic temperature, and then at a lower character- 
istic temperature to a super- Arrhenius (fragile) behavior. We address these transitions by studying 
the statistics of free energy barriers for different systems at different temperatures and space di- 
mensions. We present a clear evidence for changes in the dynamical behavior at the transition to 
Arrhenius and then to a super-Arrhenius behavior. A simple model is presented, based on the idea 
of competition between single-particle and cooperative dynamics. We argue that Arrhenius behavior 
can take place as long as there is enough free volume for the completion of a simple Tl relaxation 
process. Once free volume is absent one needs a cooperative mechanism to 'collect' enough free 
volume. We show that this model captures all the qualitative behavior observed in simulations 
throughout the considered temperature range. 



INTRODUCTION 



A full understanding of the complex dynamical sce- 
narios accompanying the glass transition requires a res- 
olution of the statistics of free energy barriers at any 
given temperature. Such information is very hard to 
come by. Trying to measure such free energy barriers 
using the energy landscape of a typical glass former is 
a daunting procedure that has never been achieved. In 
this paper we propose a method that provides us with a 
decent approximation of the statistics of the free energy 
barriers. We do this by measuring "quenched" rather 
than "annealed" relaxation times. The usual procedure 
for extracting a relaxation time from numerical simula- 
tions of any glassy system is the "annealed" procedure. 
One takes many realizations of a super-cooled liquid at 
a given temperature, measures a typical relaxation func- 
tion for each realization, averages the relaxation func- 
tions over the ensemble and finally extracts a relaxation 
time T a (T) from this average relaxation function by, say, 
determining when the average function reaches 1 /e of its 
initial value |il!]. The subscript a designates the 'main' 
slow relaxation mechanism in glassy dynamics, and is the 
only one that is typically resolved in numerical simula- 
tions. This procedure is automatically performed in ex- 
periments where the macroscopic nature of the relaxing 
systems provides self-averaging resulting in 'smooth' re- 
laxation functions that hardly fluctuate from realization 
to realization. In this paper we propose that significant 
simulational insight on glassy relaxation can be gained 
by adopting a "quenched" procedure. In this procedure 
a relaxation time r a is extracted from each and every re- 
alization by determining when the appropriate relaxation 
function reaches a value of 1/e, cf. Fig. Q] upper panel. 
Finally an average is taken over the ensemble to provide 
(r a )(T). A typical such quenched relaxation time for a 
Lennard-Jones binary glass is shown in the lower panel 
of Fig. [1] Here we measured the self part of the interme- 
diate scattering function for each realization separately 



(upper panel), i.e. the function 
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where the index i runs over half the particles in the bi- 
nary mixture (those with longer interaction length; we 
assume the mass of all the particles to be the same, 
m = 1). One observes the usual dramatic slowing down, 
such that (r a )(T) grows rapidly when T decreases, first 
in a pre- Arrhenius form (see inset), then in an Arrhenius 
form, linear in l/T, and later in a faster, super-Arrhenius 
form, which is referred to as 'fragile' behavior in the glass 
community 0. In this paper we will be interested pre- 
dominantly in the transitions to and from the Arrehenius 
regime, for which the quenched procedure is particularly 
illuminating. 

The advantage of the quenched over the annealed pro- 
cedure is that it allows a particularly transparent treat- 
ment of the statistics of free energy barriers. Having 
an ensemble of r Q values for a given ensemble of super- 
cooled systems, we can define a typical free energy scale 
by writing for each r a value an equation 



r a = t exp (F/T) 



(2) 



where we choose units such that Bolzmann's constant 
equals unity. Inverting Eq. @, every value of r a yields 
a value of F. The free energy has in principle a contri- 
bution from an energy barrier E and an entropic contri- 
bution stemming from a degeneracy g, 



F = E-T\ng 



(3) 



Of course, this definition makes sense when the dynam- 
ics involves escaping an energy barrier. The pref actor to 
is by definition the relaxation time inside the confined 
state at every temperature T, and we assume that it is 
fixed for all the realizations in the ensemble: 



r = c/VT . 



(4) 
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FIG. 1: Color online. Upper panel: the time dependence 
of ten realizations of the unaveraged intermediate scattering 
function Eq. JT]), f° r a range of temperatures (decreasing from 
left to right). By definition the relaxation time r a is the time 
at which the function reaches the value of 1/e. Notice the 
increase in the dispersion in T a when temperature is lowered. 
Lower panel: the quenched average relaxation time (r) in a 
log-lin plot vs 1/T, compared to an Arrhenius temperature 
dependence. Inset: the pre- Arrhenius regime at high temper- 
atures. The simulation for this figure were done using the bi- 
nary Lennard- Jones model whose inter-particle potential can 
be found for example in 



The dependence on 1/%/T stems simply from the esti- 
mating v, the typical particle velocity, from equipartition 
as v 2 — T and tq is then a typical length-scale (a cage 
magnitude) over v. For our purpose the constant C is 
determined by requiring that the plot of t q /to — > 1 as 
T — > oo, see Fig. [TJ(red cross) and is 0(1). 
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FIG. 2: Color online: Distribution functions P(F; T) in the 
Arrhenius regime of the Lennard Jones system in dimension 
d — 2, shown for a range of temperatures. Neither the peak 
nor the width of the approximately Gaussian distribution de- 
pends on temperature in this regime, meaning for this relax- 
ation process g = 1. 




FIG. 3: Color online: Distributions P{F; T) for the purely 
repulsive system (A), the Lennard- Jones model in d = 3 (B), 
the Kob- Andersen model (C) and the polydisperse model (D) 
in the Arrhenius regime. The distributions P(F; T) for differ- 
ent temperatures overlap throughout the Arrhenius regime. 



II. THE ARRHENIUS REGIME 

Our first observation is that in the Arrhenius regime 
the distribution of free energies F is always Gaussian, 
and is always very sharply peaked. Since F in the Ar- 
rhenius regime is temperature independent, we conclude 
that g = 1, and F = E. As an example we show again 
the binary Lennard- Jones system for which we have used 
1000 realizations for each value of T to generate the dis- 
tributions of F shown in Fig. [2] We conclude that in the 
Arrhenius regime to a very good approximation 
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with the peak value (F) and the width a unchanged 
throughout the Arrhenius regime. In the present regime 
we denote (F) = E a . It is noteworthy that the distri- 
butions are also very sharply peaked; to see this note 
that the quenched average can be now computed from 
the Gaussian integral as 
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For the data shown in Fig. g\E a = 0(1) and a ss 0.05. 
Thus for the Arrhenius range of temperatures T = 0(1) 
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the contribution of the width a 2 ps 2.5 x 10~ 3 is of the 
order of 0.1% in Eq. Obviously, for larger systems we 
expect a to be even smaller, tending to zero for N — > oo. 
It appears that all the relaxation events contributing to 
the Arrhenius regime belong to a tight group of similar 
events with a very well defined energy barrier E a and a 
negligible dispersion. In fact, we find that this conclusion 
is not particular either to the Lennard- Jones system or 
to d = 2. In Fig. [3] we show similar results for other glass 
formers with A) purely repulsive potential Q at d = 2, 
B) the Lennard-Jones system Q at d = 3, C) the Kob- 
Andersen model [5[ at d = 2 and D) the polydisperse 
model [6] at d — 2. The physical interpretation of these 
observations is discussed below, cf. Sect. I Villi 




III. THE PRE- ARRHENIUS REGIME 

At high temperatures we observe (see for example Fig. 
U) that P(F) does not have a fixed dispersion (it widens 
when T increases) and its peak moves to lower free en- 
ergy values. This observation is generic for all the models 
studied here. The physical interpretation of this obser- 
vation is again discussed in Sect. IVIIII 



IV. TRANSITION TO SUPER- ARRHENIUS 

Obviously, one might expect that when T reaches a 
value of T ps Const x E a something interesting should 
happen. The physical relaxation mechanism whose en- 
ergy barrier is E a should run out of steam, and other 
mechanisms, if they exist, should emerge. This phe- 
nomenon is beautifully seen in the Lennard-Jones system, 
see Fig. [4] Indeed, at 1/T ps 1.5 the Arrhenius process 
begins to disappear and we observe at first a gradual 
increase in both (F) (the peak position) and its width 
a. In the system for which N = 6400, we find that at 
1/T ss 2.7 the high temperature Gaussian peak disap- 
pears in favor of a second Gaussian peak which appears 
to replace it. We will denote by T* the temperature at 
which both process coexist with equal probability. The 
Gaussian of the second process marches to the right with 
increasing both of (F) and a. The transition region is 
depicted in Fig. 2] and is shown in more detail in Fig. [5] 
For the systems of the size studied in these simulations, 
by the time that the second relaxation mechanism takes 
over, the contribution of the variance a to the quenched 
average relaxation time Eq. ^ is no longer negligible. 
Thus for our example at T = 1/3 (F) « 2 and a ps 0.25 
such that a 2 / ' (F)T m 0.1. At lower temperature the dis- 
persion effect will contribute more and more to the relax- 
ation time. Of course, also here we expect a to decrease 
to zero when N —> oo. It is useful however to keep the a 
dependence for the simulations at hand. 

We should note that Eq. ([7]) remains valid throughout 
the temperature range. We can rewrite it, denoting Xi = 
(Fi)/T and introducing ua for the weighted sum of the 



FIG. 4: Color online: Distributions P(F; T) for the Lennard- 
Jones system at d = 2 with N = 6400 for the whole available 
temperature range. 
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FIG. 5: Color online: The transition in relaxation mechanism 
shown in Fig. [3] as seen by the replacement of one Gaussian 
distribution P(F) by another. 



two different mechanisms in the form 
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We test this formula for the data corresponding to the 
distributions found in Fig. [4] and demonstrate the per- 
fect agreement in Fig. [51 Note that the relaxation time 
varies in this range of temperatures by close to 5 orders 
of magnitude. 
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FIG. 6: Color online: A test of Eq.©. Shown are data for the 
binary Lennard- Jones model at d=2 and at different densities. 
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FIG. 7: Color online. Upper panel: Usual plots of 
the non-dimensional relaxation time in logarithmic scale 
vs. 1/T for five different values of the density, p = 
0.80,0.85,0.90,0.95, 1.00. Lower panel: The very same data 
collapsed by non-dimensionalizing the abscissa by the energy 
scale ep x . 



V. SCALING 

In Ref. 0] it was shown that relaxation times of models 
of the type discussed above exhibit density scaling. For 
different densities p plots of ln(r a )/r vs. 1/T fall on dif- 
ferent curves [HQ. Once plotted as a function of ep x /T 
with an appropriate value of x the data collapses onto 
one curve. Density scaling is easiest to justify when one 
uses a purely repulsive potential of the form 4>{r) ~ er~ v 
where e sets the energy scale of the interparticle poten- 
tial and r is measured in dimensionless units. Then the 
only relevant length scale £ in the system is determined 
by the density as £ ~ p" x l d where p is the density in 
dimensionless units. Obviously, when we plot ln(r)/ro 
vs. 1/T (see upper panel of Fig. [7]) we present a di- 



FIG. 8: Color online: Demonstration of the linear dependence 
of the mean Arrhenius energy E a on the depth of the potential 

£. 



mensionless number in terms of a dimensional quantity 
which is inappropriate. To correct this, we need to non- 
dimensionalize 1 /T by multiplying it by an energy scale. 
The only energy available is <j>(r), and we write 



(r) ~ ef 



ep 



,y/d 



(9) 



Thus by presenting ln(r)/ro vs. ep v ^ d /T we should 
find data collapse as discussed in Ref. p}. It is well 
known that the argument remains approximately valid 
also when the potential is not a pure power law but 
with an exponent x that needs to be found by collapsing 
the data [1(| • In the second panel of Fig. [7] we show 
the data collapse when plotted properly, with the den- 
sity varying in the range [0.80,1.00]. Additionally, in 
these systems the inter-particle potential possesses only 
one energy scale, which is nothing but the depth of the 
inter-particle potential, denoted as e. This energy scale 
must determine the typical barrier height that needs to 
be surmounted in the Arrhenius regime. Accordingly, we 
expect E a to be proportional to e. This expectation is 
fully supported by the data shown in Fig. [8]where for two 
models (binary Lennard- Jones in 2 and 3 dimensions) we 
have changed e and measured E a . 

Our additional contribution to the discussion of density 
scaling is in pointing out that it is obeyed separately 
by the two contributions to Eq. (|7|). In the upper row 
of Fig. IH1 we show that (F) exhibits a very nice data 
collapse when data for different densities are re-plotted 
as required. In the lower row we show the same for a. 
We will use these scaling functions below in providing a 
model for the observed glassy dynamics. 



VI. SYSTEM SIZE DEPENDENCE OF T* 

An obvious worry about our method of analysis is that 
the dispersion a depends on the system size, and there- 
fore the temperature T* at which we observe the transi- 
tion to fragility may depend on the system size. This is of 
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FIG. 9: Color Online: Data collapse for the normalized energy 
and dispersion for the same densities as in Fig. [7] 
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FIG. 10: Color online: Upper panel: values of T* for different 
system sizes. Lower panel: the convergence rate of T* to the 
its asymptotic value T** — 0.27. 



VII. COOPERATIVE MECHANISM 

The purpose of this section is to establish one of the 
main claims of this paper, i.e. that with the transition to 
the super- Arrhenius behavior one observes a change from 
a local to a cooperative relaxation mechanism. At tem- 
peratures below T* one expects a cooperative relaxation 
process to play an important role; the proposition of this 
paper is that this is mainly due to the disappearance 
of the free volume that was available in the Arrhenius 
regime (see quantitative model below), requiring now a 
cooperative motion of many particles to 'collect' enough 
free volume to allow a relaxation step. Although such a 
picture is current among some practitioners in the glass 
community, to the best of our knowledge it had never 
been shown explicitly. We use in this section scaling ar- 
guments to establish this picture convincingly. 

Adding the system size dependence to our typical free 
energy, we write now 



F(N,T)=E(T)-T\og[g(N,T)}. 



(11) 



We expect E(T) to be determined by the energy land- 
scape and not to depend strongly on the system size. En- 
ergy barriers are mostly sensitive to local arrangements 
of particles which are not system size dependent. On 
the other hand \og[g(N,T)] is the entropic contribution 
to the free energy barrier. Our data showed that in the 
Arrhenius regime the degeneracy factor g was of the or- 
der of unity; here we will show that our data strongly 
supports a cooperative process in the super- Arrhenius 
regime, mainly due to the system size dependence of the 
degeneracy factor. 

The flip side of having a cooperative mechanism is that 
there should exist a typical length £(T) which measures 
the degree of cooperativity and is increasing when the 
temperature decreases. Consider then a system of N 
particles in the super-Arrhenius regime, associated with 
a typical scale £,{T), and contained in a cubic box of size 
L at some temperature T . As long as the system is small, 
i.e. L <§; £(T), the whole system needs to cooperate in 
order to relax. Accordingly we expect the degeneracy 
factor to grow extensively like g(N,T) ~ N. When the 
system is large enough, i.e. when L ^> £ (T) we expect to 
find g(N,T) ~ £ d , not changing with the system size. In 
other words, 



course correct. Nevertheless we discover that the value of 
T* converges to a value T** in the thermodynamic limit. 
The data is shown for the binary Lennard- Jones model 
in Fig. 1101 From the data in the lower panel one can 
conclude that for this model 



T* - T* 



0.25 



(10) 



We did not perform an exhaustive analysis of the expo- 
nents of this law of convergence for the different models. 



](N,T) 



({)) d , for L«£(T) 



for L>£(T), 



(12) 



where a is typical inter particle distance. Accordingly the 
free energy should assume the form of a scaling function 
of L/£(T) as 



F(N,T) = E(T)-T log 
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FIG. 11: (Color online) Upper Panel: The position of the sec- 
ond peak in P(F, T) is plotted as a function of system size for 
different temperatures in the range [2.52,3.00]. Lower Panel: 
Scaling plot (see text for details) to highlight the logarithmic 
dependence on system size. 

where A is a temperature independent proportionality 
factor. The scaling function f(x) has the following 
asymptotic dependence 

Considering the small system limit, we conclude that 
there the free energy barrier should exhibit a logarith- 
mic dependence on the system size. This expectation is 
supported by the data shown in Fig 111! In the upper 
panel we show the system size dependence of the free en- 
ergy value where the distribution function P(F, T) has a 
second peak, for different temperatures. With increasing 
system size these values tend to saturate to an asymp- 
totic value. Note that the dependence on the system size 
seems to be decreasing with increasing temperature, in- 
dicating the existence of a cooperative length scale. We 
can not go to high temperatures in this study because 
at higher temperature the extraction of the second peak 
position become difficult as its height goes down drasti- 
cally, merging with the Arrhenius peak. We focus on data 
where the position of the second peak can be estimated 
accurately. To show the logarithmic dependence on the 
system size we collapsed the data together by rescaling 
the free energy axis by £ 5 with § ~ 1.2 and the x axis by 
£(T) d . For simplicity we fitted here a guess function for 
}(T) in the form f(T) - 1/T 18 [11]. In the lower panel 
of Fig [TT] one can see the convincing logarithmic depen- 
dence of the free energy on the system size for different 
temperatures. 



VIII. INTERPRETATION AND DISCUSSION 

The data shown above indicates two changes in the 
relaxation mechanism between the pre- Arrhenius and the 
Arrhenius regime and then between the Arrhenius and 
the fragile regime. We propose that the pre-Arrhenius 
regime is a fluid regime, in which a well-defined cage 
is only beginning to form around every particle. It is 
relatively easy to break this ill-formed cage, and there 
are many ways to do it. Accordingly the free energy 
needed is made of a low value of E and a finite value of 
g such that when T is increasing the Gaussian peak that 
we observe marches to the left and becomes broader. 

In the Arrhenius regime a well defined cage has formed 
around each particle, and this cage is a more or less reg- 
ular cluster of particles aggregating around each center 
particle. We expect a roughly constant energy E to suf- 
fice to break this kind of cage and the width of the dis- 
tribution of F should be relatively narrow as discussed 
below. Towards the end of the Arrhenius regime a second 
and maybe third layer of next-nearest neighbors begins to 
form around a center particle, and therefore it becomes 
more and more difficult for the single-particle relaxation 
to take place. The energy barrier increases but because 
of the fluidity of the second and third layer also the width 
of the distribution starts to increase. 

Eventually a new, cooperative mechanism must be fa- 
vored for dynamical relaxation. Note that the appar- 
ent mean energy of activation (F) increases when the 
temperature decreases below the Arrhenius regime. This 
must mean that the single particle relaxation mechanism 
that is operative at the Arrhenius regime is no longer 
available at lower temperatures, since if it were available 
it would have been selected. In the temperature range 
where both mechanisms coexists we indeed find the asso- 
ciated (F) and a to be similar. We thus offer a model of 
the observed phenomenology on the the following basis: 

(i) In the Arrhenius regime the relaxation process is an 
elementary event like a T± process which is characterized 
by a relatively sharp free energy cost F\, 

F 1 = E a -T]ng ll (15) 

where g± is the number of available inequivalent Tl pro- 
cesses. This energy cost E a is of the order of breaking 
through a local cage and the process involves only a few 
particles. Due to the amorphous nature of the super- 
cooled liquid E a is a stochastic variable but the dispersion 
around E a is relatively small. We denoted the variance 
in the distribution of E a as a 2 and noted that in this 
regime a/E a <C 1. 

(ii) The fact that Fi does not change with the temper- 
ature in the Arrhenius regime indicates that g± is of the 
order of unity. 

(iii) The precise reason for the blocking of this elemen- 
tary relaxation channel when temperature is lowered is 
not known. We will propose, however, that this block- 
ing is due to the disappearance of available free volume 
per particle v/(T,P) [12j. In other words, denoting by 
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v m the minimal average volume that every particle occu- 
pies, there must be additional free volume Vf for the T\ 
process to occur. 

Of course, the average volume per particle, which is 
the inverse of the density, v(T,P), is determined by the 
equation of state which varies from one material to the 
other. We will define Tf as the temperature where the 
free volume disappears, and then v(Tf, P) = v m . For the 
purposes of the present model we estimate the average 
volume per particle from the first order Taylor expansion: 



v(T,P)&v(T f ,P) + 



dv 
df 



(T — T- 



f) 



(16) 



where the partial derivative is computed at T = Tf. Ac- 
cordingly, since v(Tf,P) — v m , 

Vf(T,P) « (jjf) ( T ~ T f) .provided v > v m (17) 



Vf(T,P) m otherwise 



(18) 



The probability of having a local free volume Vf is 
estimated from maximum entropy considerations as 



P (vf) = — ex P -(vf/vf) . 



(19) 



Consequently the probability that vt > v m , which is the 
probability that a Tl process could occur, is 
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FIG. 12: Color online. The model prediction for the develop- 
ment of the statistics of F throughout the temperature range, 
from pre- to post-Arrhenius behavior. The parameters were 
chosen to fit the the data of the Lennard- Jones system in 2d, 
and accordingly the results are in excellent agreement with 
the numerical results shown in Fig. f4] For much larger sys- 
tem sizes we expect the dispersion a to tend to zero, but this 
affects neither our model nor the conclusions. 



because the T\ process is no longer available or it requires 
cooperation from the surroundings to be realized. Also 
an increase in the number of particles involved must in- 
crease 52 simply due to the number of ways that these 
particles can be chosen. Thus we write in the fragile 
regime 



P(vf > iv, 



P(vf)dv f 



(20) 



Finally, we write Vf = v — v m = av m (T — Tf), where a is 
the isobaric coefficient of thermal expansion at T = Tf , 



1 



With this we can write 



dv 
~dT 



(T = T f ) 



P(vf>v m )=cxp [-cTf/(T-T f )} 



(21) 



(22) 



where the dimensionless constant c = 1/aTf. Note that 
we can immediately estimate the probability of seeing the 
cooperative mechanism (which we will refer to as the T^ 
process) as simply 1 — exp [—cTf/(T — Tf)]. Therefore we 
can estimate T* where the two processes are equiprobable 
from 



exp [-dT f /{T*-T f )]n 1/2 



(23) 



and therefore T* w (1 + c/ln2)X/. 

(iv) In the fragile regime we propose that Eq. (fT5)l 
retains its form, but now with a temperature dependent 
E 2 (T) and a sizeable degeneracy gi (T). Our picture is 
that of a relaxation process that involves more and more 
cooperative motions of a larger and larger number of par- 
ticles. This obviously results in an increasing energy bar- 
rier Ei{T) which is indeed higher than E a but selected 



F 2 {T) = E 2 {T)-T\ng 2 {T) 



(24) 



We now consider the emergence of the T 2 process as a 
consequence of the disappearing of the T\ process as is 
presented for example in Fig. [4] We write 



P(F.T) 

Z(T) 



W 1 (F,T) + W 2 (F,T) 



(25) 



Z(T) 

J \W\(F, T) + W 2 (F,T)}dF . (26) 



Here W\ and W 2 are the weights of the two relaxation 
channels respectively. In the Arrhenius regime W\ is 
composed of three factors; remembering that g\ ~ 1, 
we write in that regime 



W l (F,T)=e 



- J-cT f /(T-T f )] 



1 



V27 



x e 



-F/T 



(27) 

The first factor is the probability to have a T\ process, 
Eq. [521 The second factor is the probability to select a 
free energy barrier of magnitude F. The third factor is 
the probability to overcome that particular barrier. In 
the larger temperature regime that includes the pre- and 
post-Arrhenius regime we generalize Eq. (l2"T|) to read 



Wx{F,T) 



o [-cT f /(T-T f )], 



1 



- [F -(F 1 >(T)]^ 



xe 



-F/T 



(28) 
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where of(T) « er 2 + K(T - T AlI ) 2 . This fit is supported 
by the scaling function shown in Fig. [5] in the left lower 
panel, and T AlI is the temperature of the minimum of 
scaling function. 
Similarly, 



W 2 (F,T) 



U - e [- cT //( T - T /)]| -[ f -<f2>(T); 



x e 



-F/T 



(29) 

In Fig. Q2] we present the prediction of Eq. (|2"o]) as a 
function of temperature. We chose the various param- 
eters to agree with the functions shown in Fig. [9] with 
c = 0.6 and T* = 0.55. The model indeed reproduces 
very well the main observations discussed above. We 
see the pre-Arrhenius behavior and its transition to Ar- 
rhenius behavior where expected. We also observe the 
rapid decline of the Arrehnius gaussian with the concur- 
rent increase of the fragile gaussian which is marching to 
the right with increasing dispersion. We should stress 
that this simple model is not meant to represent the 



full complexity of the observed processes, especially at 
low temperatures where the cooperative dynamics reigns 
supreme; the point of this simple model is to stress the 
fundamental observation of this paper, i.e. that a com- 
pact and well defined relaxation mechanism appears to 
operate in the Arrhenius regime, whereas this simple 
mechanism is blocked once temperature reduces, giving 
rise to another mechanism which is presumably coop- 
erative and much more involved. For all that we know 
one may encounter additional relaxation processes at still 
lower temperatures which are currently not available to 
numerical simulations. 
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